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It has been found in recent heavy-ion experiments that the second and the third flow harmonics 
of direct photons are larger than most theoretical predictions. In this study, I construct effective 
parton phase-space distributions with in-medium interaction using quasi-particle models so that 
they are consistent with a lattice QCD equation of state. Then I investigate their effects on thermal 
photons using a hydrodynamic model. Numerical results indicate that elliptic flow and transverse 
momentum spectra are modified by the corrections to Fermi-Dirac and Bose-Einstein distributions. 
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I. INTRODUCTION 

The discovery of collective dynamics of the QCD media 
created in BNL Relativistic Heavy Ion Collider (RHIC) 
and CERN Large Hadron Collider (LHC) implies early 
equilibration of low-momentum gluons and light quarks, 
and realization of a quark-gluon plasma (QGP) fluid 
[1, 2]. Dynamical description of a QCD system in the 
vicinity of the quark-hadron crossover is generally a diffi¬ 
cult issue. The fluidity of the system, on the other hand, 
is believed to allow an effective description of the bulk 
dynamics of heavy-ion collisions [3] where aside from the 
initial conditions the information of QCD comes into the 
formalism through the equation of state (EoS), the trans¬ 
port coefficients and the chemical reaction rate. Particle 
spectra and flow harmonics [4, 5] of hadrons are reason¬ 
ably well-explained in the framework of hydrodynamics 
with small viscosity. Direct photons, on the other hand, 
are found to be non-trivial, because elliptic flow V 2 [6, 7] 
and triangular flow [8] of direct photons undershoot 
those measured in the experiments by a factor when es¬ 
timated in conventional models. Here direct photons are 
defined as the photons which does not originate from 
hadronic decay. Initial hard photons, named prompt 
photons, do not have azimuthal anisotropy, but medium- 
induced soft photons, known as thermal photons, can 
inherit anisotropy from the medium. The apparent dis¬ 
crepancy poses theoretical challenges to the community 
[9-31]. ^ 

Despite the experimental evidences that the system 
is likely to be strongly-coupled, heavy-ion models often 
employ the phase-space distributions of ideal gases. In 
general the equilibrium phase-space distributions in an 
interacting system should be subject to non-ideal correc¬ 
tions. They are known to play important roles in phe¬ 
nomenology in non-relativistic systems - carbon dioxide, 
for example, has a solid phase due to van der Waals force 
which arises from interaction corrections. It should be 
noted the deviations from the ideal gas picture is a differ¬ 
ent concept than viscous corrections because the system 
is in local equilibrium with no entropy production. 
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In heavy-ion collisions, hydrodynamic flow is converted 
into hadrons at freeze-out using the one particle distribu¬ 
tion so that energy-momentum tensor is conserved on the 
hypersurface using relativistic kinetic theory [32]. This 
is supported by the fact that the lattice QCD EoS agrees 
with that of hadron resonance gas well - but only ap¬ 
proximately - below the crossover temperature. Also 
the second-order transport coefficients which appear in 
causal viscous hydrodynamic equations [33] are often an¬ 
alytically derived in kinetic theory and the energy density 
and pressure in the expressions are substituted by the hy¬ 
drodynamic ones, which partially introduces the effects 
of interaction to the transport coefficients. 

For photons, the parton distribution functions are used 
in thermal photon emission rates. Oftentimes calcula¬ 
tions are performed simply assuming Fermi-Dirac and 
Bose-Einstein distributions for quarks and gluons, re¬ 
spectively, in the QGP phase, though lattice QCD re¬ 
sults seem to deviate from the results implied from those 
assumptions even at relatively high temperatures [34]. 
Considering the recent status of the photon V 2 (and V 3 ) 
puzzle, the real gas effects should be worth-investigating 
because the reduction of the effective degrees of free¬ 
dom at high temperatures could suppress early-time pho¬ 
tons with small anisotropy while hadronic photons with 
larger anisotropy would not be affected much due to the 
aforementioned agreement of the kinetic theoretical es¬ 
timations and the lattice QCD data. There have also 
been important studies to improve thermal photon mod¬ 
els through the phase-space distribution. A few recent 
examples include viscous corrections to the distribution 
[35, 36] and semi-QGP effects [30]. 

In this paper, I discuss the effects of the interaction in 
parton phase-space distributions on heavy-ion photons 
based on the quasi-particle parametrization. There are 
many variants of the quasi-particle models [37-43] . Here 
the model is phenomenologically constructed so that it 
reproduces thermodynamic variables of the lattice QCD 
equation of state. Then I perform numerical estima¬ 
tions of thermal photon elliptic flow and particle spectra 
by consistently using the same EoS in a hydrodynamic 
model and show that modification of the emission rate 
of the QGP photons leads to visible corrections to the 
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observables. 

In Sec. II, the quasi-particle model is developed based 
on the lattice QCD results. The model for the estima¬ 
tions of thermal photons and background hydrodynamic 
flow is presented in Sec. III. Sec. IV is devoted for numer¬ 
ical results. Sec. V presents conclusions and discussion. 
The natural units c = h = ks = ^ and the Minkowski 
metric = diag(-|-, —, —, —) are used in the paper. 


II. QUASI-PARTICLE MODEL 


I employ the quasi-particle picture to estimate the ef¬ 
fects of in-medium interaction on thermal photon V 2 ^ be¬ 
cause the photon emission rate is usually given as a func¬ 
tional of distribution functions. Here the quasi-particle 
model is constructed so that it is compatible with lattice 
QCD results. The effects of interaction is represented 
by the (self-)interaction term Wlff{p,T) where p is the 
four-momentum and T is the temperature. The effective 
one-particle distribution is 

= exp(w,/r)±r 


and the grand-canonical partition function Zi in a loga¬ 
rithmic form is 


In Zi = ± V 


9idp^ 

(27r)3 


In 


1 ± exp 



V 
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where V is the volume, gi is the degeneracy, uji = 
yjp^ mf + IVgg is the effective energy density, and i 
is the index for particle species, i.e., quarks and gluons. 
Here the numbers of quarks and antiquarks are assumed 
to be the same, i.e., the net baryon current is not con¬ 
sidered. I consider the light quarks u, d, and s for the 
equilibrated quark components. The sign of quantum 
statistics is -I- for fermions and — for bosons. <I>(T) is 
the background field contribution, which is determined 
by the thermodynamic consistency condition [38] 

^1 /O', 

dT J (27r)3 dT ^ ’ 

and the condition P(T = 0) = 0. This correspond to the 
bag constant in the bag model [44] but here it is temper¬ 
ature dependent. Wl^ is assumed to include the infor¬ 
mation on interaction and self-energy correction to the 
hamiltonian and can be regarded as an effective chemi¬ 
cal potential as well as thermal correction to the particle 
mass. Energy density and hydrostatic pressure are writ¬ 
ten as 




FIG. 1. (Color online) Lattice data fit of the effective inter¬ 
action Weff as a function of the temperature. Dashed line 
denotes the region below the crossover. 
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where /3 = 1/T and 4> = 4*^. The last line of Eq. (4) is 

obtained by using the consistency condition (3) and that 
of Eq. (5) by integration by parts. It should be noted that 
the trace anomaly = e — 3P is no longer vanishing 
in the massless limit due to the presence of the effective 
interaction. The entropy density is given through the 
thermodynamic relation s = dP/dT = {e + P)/T. 

Determination of the effective interaction contribution 
W^ff is generally a non-trivial issue. Here it is constrained 
with the (2-|-l)-flavor lattice QCD equation of state [45] 
assuming it depends only on the temperature. For sim¬ 
plicity and the lack of additional constraints I assume 
that Weff is common for all partons in this paper to see its 
qualitative effects. It should be noted that a thermody¬ 
namically consistent formalism yields the energy density, 
the pressure, the entropy, and the trace anomaly simul¬ 
taneously. One can use the entropy density to determine 
Weff as the quantity is independent of 4’, and then deter¬ 
mine the background contribution by the relation (3). 

Figure 1 shows the temperature dependence of Weff 
fit to the lattice data. It has a peak structure slightly 
below Tc because one has to mimic the crossover behav¬ 
ior which does not exist in the free parton gas picture. 
One can see that the effective fugacity would be non- 
negligible because it is comparable to the temperature. 
The corresponding background contribution 4> in a di¬ 
mensionless form is shown in Fig. 2. Compared with the 
energy density and the pressure, the effective background 
contribution is small. Functional fits to Weff and $ are 
presented in Appendix. A. It should be noted that the 
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FIG. 2. (Color online) The background contribution in a 
dimensionless form as a function of the temperature. 

Dashed line denotes the region below the crossover. 


parton picture is at most valid above the crossover tem¬ 
perature. 

Although quasi-particle interpretations are commonly 
employed in hadron physics, I emphasize that the pur¬ 
pose of considering it is to improve the ideal gas distri¬ 
butions often assumed in thermal photon emission rates. 
The real QCD system is not necessarily a relativistic 
gas of quasi-particles especially near the crossover. One 
can also construct similar effective distributions for the 
hadronic phase but it is not considered here because the 
difference between the EoS of hadronic resonance gas and 
that of lattice QCD is indicated to be relatively small as 
mentioned earlier. 


III. THE MODEL 

I develop the numerical model to estimate thermal 
photon elliptic flow and transverse momentum spectra 
with and without the effects of interaction corrections us¬ 
ing the relativistic hydrodynamic model and the thermal 
photon emission rates in which the effective distributions 
obtained in Sec. II are embedded. 


A. Thermal photons 


The thermal photon emission rate is affected by the 
modification of the phase-space distribution. For the 
demonstrative purpose, the hard photon emissions in 
Compton scattering and pair annihilation are consid¬ 
ered for the QGP phase. I employ the emission rate in 
Ref. [46] at the leading order in the fugacity expansion 
by substituting the parton distributions with the effective 


distributions. The emission rate reads, 


di? _ Soos 2^-B/T 

~ 97r2 


AgAg 


log 


+ A 
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Here Xq = \g = are the effective fugacities, 7 is 

Euler’s constant, and = g^T^/6 is the infra¬ 

red cut-off. 

The hadronic photons are assumed to be unaffected. 
The emission rate in the hadronic phase is employed from 
Ref. [47]. It should be noted that the photon emission 
rate near the crossover can be in principle non-trivial 
[28] . Here they are simply interpolated with a hyperbolic 
function as 




where c(T) = {1 -|- tanh[(T — Tc)/AT]}/2 with the con¬ 
necting temperature Tc = 0.17 GeV and the crossover 
width AT = 0.017 GeV. 


B. Hydrodynamic flow 


The flow and the temperature profiles are calculated 
using the ( 2 -|-l)-dimensional ideal hydrodynamic model 
assuming boost-invariance in the space-time rapidity di¬ 
rection [25]. The net baryon number is assumed to be 
vanishing. The hydrodynamic equation of motion is 
energy-momentum conservation d^T^'^ = 0 where 

T^'^ = {e + P)u>^u'' - Pg^’', ( 8 ) 


remains the same for ideal and effective distributions, 
is the flow. The medium property is determined by the 
EoS, which provides an additional equation for uniquely 
determining all the thermodynamic variables. Here it is 
employed from the aforementioned lattice QCD results 
[45], which of course is consistent with the quasi-particle 
description we use for the estimations of thermal pho¬ 
tons. The Monte-Carlo version of Glauber model [48] is 
used to construct a smoothed initial condition. The nu¬ 
merical code is from Ref. [27]. Here the initial condition 
is obtained by averaging over events. 

Thermal photon particle spectra can be calculated by 
integrating the photon emission rate over space-time 

dN'^ _ 7 ^4 dR'^ 

dcppprdprdy J d(l)pPTdpTdy ’ 

and elliptic flow by taking its second-order Fourier har¬ 
monics in the azimuthal angle relative to the reaction 
plane angle 41 in momentum space as 


vl{PT,y) 


fo" cos(2(^p - 4') 

dfh _ iAJ. _ 

JO d(l)pPTdpTdy 


( 10 ) 


Here pt is the transverse momentum, y is the rapidity, 
and (j)p is the angle in momentum space. The local photon 
energy is estimated as E = p^Ug. 
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IV. NUMERICAL RESULTS 

For demonstrative purposes, Au-Au collisions at 
•^/snn = 200 GeV with the impact parameter 6 = 6 fm 
are considered. The energy density is normalized so that 
the maximum energy density is e„iax = 30 GeV in the 
most central collisions at the initial time of hydrody¬ 
namic evolution, which is chosen as tq = 0.4 fm/c. The 
contribution of the photon emission above the hadronic 
freeze-out temperature Tf = 0.13 GeV is taken into ac¬ 
count. 


A. Thermal photon V2 

The differential thermal photon elliptic flow with the 
ideal and the effective distributions at mid-rapidity are 
shown in Fig. 3 for the transverse momentum window 
0.2 < pt < 5 GeV. One can see that the effective in¬ 
teraction corrections enhance v^iPT)- This can be un¬ 
derstood as follows. The emission of thermal photons 
are suppressed at early times due to the fact that the 
effective degrees of freedom in the quasi-particle system 
is smaller than that in the free gas system when the sys¬ 
tem is in the QGP phase. The thermal photons emitted 
at late times have larger elliptic flow because azimuthal 
anisotropy in the medium develops along with time evo¬ 
lution driven by the pressure gradients. Thus when early- 
time contribution becomes effectively small, the overall 
elliptic flow can become large. It should be noted that 
unlike hadrons, photons are emitted at each space-time 
point and they all contribute to the final spectra because 
their interaction with the bulk medium would be weak. 

The effect of the non-ideal gas corrections is consistent 
with the recently-observed large photon vJipx), though 
the numerical estimations suggest that the effect could 
be part of the cause for the excessive photon anisotropy 
but would not be large enough to be the sole reason with 
the current hydrodynamic parameter sets. 


B. Thermal photon px spectra 

Figure 4 shows the pr spectra of thermal photons at 
mid-rapidity with ideal and effective distributions, re¬ 
spectively. The photon spectrum with the interaction 
effects is naturally suppressed by the suppression of ini¬ 
tial QGP photon emission. This implies that additional 
photon emission source might be required to explain the 
spectra when the interaction effects in the phase-space 
distributions are properly taken into account. Similar 
discussion can be found, for example, in the system with 
incomplete quark chemical equilibration [25]. The mag¬ 
nitude of the spectra suppression and V 2 enhancement 
effects are sensitive to the relative ratio of QGP photons 
to hadronic ones. 

The proper time r dependence of the photon emission 
from the medium is also investigated. Here it is defined 



FIG. 3. (Color online) Thermal photon elliptic flow with ef¬ 
fective distributions (dashed line) compared to that with ideal 
distributions (solid line). 



FIG. 4. (Color online) Thermal photon pr spectrum with 
effective distributions (dashed line) compared to that with 
ideal distributions (solid line). 


as the emission rate integrated in the spatial and the 
transverse momentum directions, 

= J dx^^^{T,x,y). ( 11 ) 

The numerical results are shown in Fig. 5 for the sys¬ 
tems with ideal and effective distributions, respectively, 
at y = 0. At the beginning photon production is reduced 
along with time evolution because the temperature de¬ 
creases rapidly due to the longitudinal expansion. The in¬ 
medium corrections suppress the emission by a factor of 
~ g- 2 Wett/T, becomes larger as it nears the pseudo- 

critical temperature. On the other hand, hadronic pho¬ 
ton contributions become larger as the system cools down 
and consequently the difference between the two systems 
disappears at later times. 
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FIG. 5. (Color online) Spatially-integrated photon emis¬ 
sion as a function of proper time with effective distributions 
(dashed line) compared to that with ideal distributions. 


V. DISCUSSION AND CONCLUSIONS 

Effects of the interaction corrections in parton distri¬ 
butions on heavy-ion thermal photons are investigated. 
This allows one to preserve the consistency between the 
hydrodynamic evolution and the photon emission rate. 
The quasi-particle model is introduced and constrained 
so that it reproduces the thermodynamic variables from 
the lattice QCD estimations. Thermal photon emission is 
estimated in numerical simulations, and it is found that 
the suppression of the QGP photons in early stages leads 
to enhancement of thermal photon V 2 - px spectrum of 
thermal photons, on the other hand, is suppressed by 
the corrections. The results imply that one has to take 
the modifications of phase-space distributions into ac¬ 
count for the quantitative description of heavy-ion pho¬ 
ton spectra and flow harmonics. Since the mechanism of 
early-time photon suppression can also lead to the en¬ 
hancement of thermal photon U 3 , it would be important 
to perform event-by-event analyses for calculating higher- 
order flow harmonics. 

One should be careful before quantitative discussion 
that prompt photons, which will reduce the overall 
anisotropy in direct photons, are not included here and 
that only hard photons above are considered for the 
QGP photons. Also though the corrections to ideal gas 
distributions would be a necessarily ingredient in the 
heavy-ion modeling and the magnitude of V 2 enhance¬ 
ment is visible in the model calculations, it would not be 
enough to fully account for the large direct photon V 2 ob¬ 
served in the experiments especially with prompt photon 
contributions. 


There would be several interesting applications of 
the quasi-particle model to heavy-ion phenomenology. 
Firstly, the formalism can provide a model equation of 
state for chemical non-equilibrated quark-gluon systems 
which is compatible with the lattice QCD EoS in the 
equilibrium limit, which is useful in studying the effects 
of quark chemical equilibration in hydrodynamic systems 
[25] . Secondly, it would be interesting to apply the effec¬ 
tive distributions for the derivation of causal dissipative 
hydrodynamic equations because some formalism, such 
as Israel-Stewart theory [33], are explicitly dependent 
on the equilibrium distribution. The transport coeffi¬ 
cients can then be fully consistent with hydrodynamic 
EoS. Thirdly, this may also allow one to formulate an 
anisotropic hydrodynamic formalism which is compati¬ 
ble with the lattice EoS by using the effective distribution 
instead of the ideal gas distribution. 
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Appendix A: Eunctional Fits for the Effective 
Distributions 

In this appendix, the polynomial fits to the in-medium 
correction parameters in the parton phase-space distri¬ 
butions are presented. The contribution of gluons and 
up, down, and strange quarks are considered. The fit to 
the effective interaction energy WeS for 0.15 < T < 0.6 
GeV is given as 

tUeff = a2t'^ + aii + ao + ^ + ^ + ^ + ^, 

(Al) 

where 02 = —1.040, oi = 2.419, oq = —1.879, a_i = 
0.8181, a _2 = -0.1796, a _3 = 0.02072, and a _4 = 
—0.0009111. Here t = T GeV/1 GeV is defined as the 
dimensionless temperature. 

Likewise, the background contribution $ would be 
given by the following functional fit: 

^ + b2t'^ + ht + bo, (A2) 

where 63 = —0.3303, &2 = 0.2424, bi = —0.04152, and 
bo = 0.001981 for the aforementioned temperature range. 
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